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A thermoacoustic compressor is capable of converting an alternating gas flow to a direct one with a large 
pumping rate on the basis of the pressure oscillation nature of thermoacoustic engines and the flow 
rectification effect of check valves. Theoretical calculations are first carried out to study the factors that 
affect the performance of the closed and open thermoacoustic compression systems. It is shown that the 
frequencies of directly connected thermoacoustic engines should avoid small integer multiple re¬ 
lationships to operate efficiently. Increasing the pressure amplitudes is beneficial for the pressure lift in a 
closed system as well as the pumping rate in an open system. A demonstrative closed thermoacoustic 
compressor was then experimentally studied. A maximum average gas pumping rate of 4.55 Nm 3 /h 
during the first 2 s of the compression process was achieved when all components were at the same 
initial mean pressure of 2.13 MPa. The maximum pressure lift reached 0.4 MPa when the initial mean 
pressure was 2.4 MPa. It was found that the pressure lifts were roughly proportional to the pressure 
amplitudes. Due to the superposition of alternating and direct gas flows, deformation of pressure 
waveforms which has a negative effect on the performance was observed. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Both alternating and direct flows are widely used in industries, 
such as in reciprocating and centrifugal compressors. In some ap¬ 
plications, the two flow types even coexist. Taking G-M cryocooler 
[1 as an example, the gas flow supplied by the scroll helium 
compressor is direct flow while what required by the cold head is 
an alternating one. The rotary valve between the helium 
compressor and the cold head converts the direct gas flow into an 
alternating one by periodically switching between the high and low 
pressure ports of the compressor. 

In recent years, the technology of flow type conversion has given 
birth to some novel promising applications. Recently, a mean flow 
acoustic engine composed of a main pipe and two single end closed 
branch tubes was proposed to induce acoustic oscillation from the 
mean flow (direct flow) in the main pipe [2-6]. Based on this 
conversion technique, aerodynamically driven thermoacoustic re¬ 
frigerators and piezoelectric generators can be developed for mean 
flow energy (e.g. wind energy, fluid energy in pipelines) exploita¬ 
tion. With the development of smart materials, piezoelectrically 
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driven hydraulic micropumps have been proposed for applications 
such as chemical process control, drug delivery system, and cooling 
of tiny electric element [7-9 . The hydraulic pump utilizes reed 
valves, diffuser nozzles and active piezoelectric discs to convert the 
reciprocating motion of a piezoelectric actuator to a pumping ac¬ 
tion. In 2004, a resonant self-pumped circulating thermoacoustic 
heat exchanger was presented by Swift and Backhaus [10] to ex¬ 
change heat with a remote heat source. Hydrodynamic asymmet¬ 
rical constrictions were adopted as an imperfect gas diode to 
generate substantial mean flow out of a thermoacoustic engine. 
Later, check valves were adopted by Gao et al. [11] to realize the 
conversion from sinusoidally alternating to direct gas flow in a non¬ 
resonant self-circulating thermoacoustic heat exchanger. 

Recently, we proposed a remarkable thermoacoustically driven 
compression effect based on the conversion of gas flow from an 
alternating state to a direct one [12-14 . Due to the unique features 
of thermoacoustic engines [15-18 , thermoacoustic compression 
systems have the advantages of simple structure, thermally driven 
nature, oil-free operation and large operating range, and may find 
applications in gas compression, power generation, and refrigera¬ 
tion systems. If an expander is introduced further, a closed cycle 
solar power generation system or cooling system can be formed 
[14]. Besides, because of the already-demonstrated feasibility of 
micro-miniaturized thermoacoustic engines 19] and check valves 
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[20-22], micro thermoacoustic compressors or pumps may find 
potential applications in MEMS field in the future. 

In the previous paper [12 , the thermoacoustic compression 
process is first proposed and a demonstrative thermoacoustic 
compressor verifies the feasibility of the concept. In order to reveal 
the operating principles and show the basic guidelines for the 
design of a thermoacoustic compression system, a detailed study of 
the working characteristics of closed and open thermoacoustic 
compression systems is carried out. The effects of operating fre¬ 
quency, pressure amplitude, phase difference, and mean pressure 
on the performance of such systems are presented in detail. 

2. Computation and analysis 

Fig. 1(a) shows the schematics of two types of thermoacoustic 
compression system, i.e., a closed and an open thermoacoustic 
compression systems. The CTCS (closed thermoacoustic compres¬ 
sion system) consists of two or more thermoacoustic engines and 
in-between check valves, as denoted in the dashed box. The num¬ 
ber of thermoacoustic engines is always one more than that of 
check valves in CTCS. The OTCS (open thermoacoustic compression 
system) consists of an arbitrary number of thermoacoustic engines 
and one more check valves. Unlike the CTCS inside which gas is 
always restricted, the OTCS continuously pumps gas from a low 


are calculated. For the OTCS, as the gas is continuously supplied to 
the LMPE and discharged from the HMPE, the mean pressures stay 
constant. Thus, the pumping rate through the check valve is what 
applications concern. Therefore, a basic OTCS unit consisting of a 
check valve and two thermoacoustic engines at constant mean 
pressures is studied, and the average pumping rate through the 
check valve is calculated. The average pumping rates in both CTCS 
and OTCS in the computations are calculated via, 
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where R g , To, and po are the gas constant, temperature, and pressure 
under standard conditions, respectively, m and r denote the total 
gas mass transported from the LMPE to the HMPE and the time 
needed to complete the compression process in CTCS. While in 
OTCS, m denotes the total gas mass transported in the calculation 
time r, which is 10 s in the calculation, m and U denote the 
instantaneous mass flow rate and volumetric pumping rate 
through the check valve respectively. 

The instantaneous mass flow rate m is determined by the 
pressure difference across the check valve Ap, and is calculated by 
Refs. [12,23]: 
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pressure gas source to a high pressure one. Check valves are set at 
the downstream and the upstream sides of the low and high 
pressure gas sources, respectively. CTCS can easily demonstrate the 
working principle of the thermoacoustic compression, while OTCS 
is more likely to be used in industry where continuous gas 
compression and pumping are required. The arrows in the figure 
show the directions of the gas flow. Due to the thermoacoustic 
conversion effect, the pressures p H and pi in the HMPE (high mean 
pressure engine) and the LMPE (low mean pressure engine) oscil¬ 
late periodically with mean pressures pn.m and pL, m . respectively, as 
shown in Fig. 1(b). As denoted by the shadows, the pressure pi in 
the LMPE can be higher than the pressure pn in the HMPE at some 
moments due to pressure oscillation. The check valve enforces a 
one direction gas flow from the LMPE to the HMPE when the 
pressure difference Ap = pi - pn exceeds the opening pressure 
difference Ap c of it. This is the main operating principle of a ther¬ 
moacoustic compressor. 

To simplify calculation, both the basic CTCS and OTCS units with 
nitrogen as working gas are studied. The basic CTCS unit consists of 
two thermoacoustic engines and a check valve. The volumes of the 
LMPE and the HMPE are 6.3 L and 5.3 L, respectively. The initial 
mean pressures of the LMPE and the HMPE are both set at 2.4 MPa 
unless otherwise stated. Due to the thermoacoustic compression 
effect, the gas in the LMPE will be transported to the HMPE until the 
transient pressure difference Ap across the check valve can not 
exceed Ap c anymore. The final pressure lift Ap m and pumping time 


where p, A, I< and F mu it denote gas density, flow area of the check 
valve, local loss coefficient and empirical multiplier for viscous 
pressure drop, respectively. R, generally a positive number «1.0, is 
the ratio of fully-closed to fully-opened flow area of the check valve. 
In the present calculations, I<, F mu i t , and R are set as 2.8,1, and 0.01, 
respectively. The opening pressure difference Ap c and the maximal 
inner diameter of the check valve are set at 0.01 MPa and 6 mm, 
respectively. To simplify the analysis, it is assumed that the pres¬ 
sure waveforms are sinusoidal and the response time of check valve 
is short enough to be ignored. The pressure waveform of the HMPE 
leads that of the LMPE by a phase 6. The time step is set at 1 x 10 6 s. 

The direct driving sources of the gas compression and pumping 
in thermoacoustic compression systems are the pressure oscilla¬ 
tions in thermoacoustic engines excited by the input thermal en¬ 
ergy. With the thermoacoustic energy conversion effect, the input 
thermal energy is converted into acoustic power, which is then 
consumed mainly in the following ways: gas pumping through 
check valves, viscous and thermal-relaxation effects in thermoa¬ 
coustic engines. Thus, the efficiency of the whole system is the 
product of thermal-to-acoustic efficiencies of thermoacoustic en¬ 
gines and the compression efficiency of the compression process. 
The thermal-to-acoustic efficiencies are strongly influenced by the 
design and the operation of the thermoacoustic engines. Typically, 
the thermal-to-acoustic efficiency of a well-designed standing- 
wave thermoacoustic engine can reach about 0.15-0.20 [24 . As 
this study mainly focuses on the characteristics of gas compression 
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Fig. 1 . Schematic for showing working principle of a thermoacoustic compressor. 


processes, the thermal-to-acoustic efficiencies of the thermoa¬ 
coustic engines are not discussed thoroughly here. For simplicity, 
the compression efficiency tj of the OTCS is defined as the theo¬ 
retical minimum work requirement W\ divided by the actual 
acoustic power consumption Wd, which is actually an indicator of 
the perfection degree of the gas compression process, 


7) = 


w { 

w 



The theoretical minimum work required is given by the 
isothermal compression process, 


W ]= mR g T 0 \n^. (4) 

PL,m 

Neglecting the acoustic power dissipation in the connecting 
tubes, the actual acoustic power consumption Wd in a thermoa¬ 
coustic compression system can be expressed as 

W d =E cv + E L + E H , (5) 

where E cv , El and Eh are the acoustic power dissipation terms in 
check valve, LMPE and HMPE. The acoustic power dissipation in 
check valve E cv can be calculated by 

0 

For standing-wave thermoacoustic engines, the acoustic 
power dissipations in thermoacoustic engines can be estimated by 
Ref. [25] 
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in which pi and Eli denote local pressure and volume flow ampli¬ 
tudes along the axial direction. f K and / v are thermal and viscous 
functions [25 . A g is the cross-sectional area of the gas flow channel. 
(D is the angular frequency of the engine. The distributions of the 
pressure and volume flow amplitudes in a standing-wave ther¬ 
moacoustic engine can be calculated via thermoacoustic transfer 
matrix method [26,27]. When the dimensions of thermoacoustic 
engines and the working conditions of a specific OTCS are given, the 
compression efficiency ij can be calculated. 


2 A. CTCS (closed thermo acoustic compression system) 

2 A A. Effect of frequency 

The effect of the operating frequency on the pressure lift Ap m of 
the CTCS is shown in Fig. 2./l and/n are the operating frequencies of 
the LMPE and the HMPE, respectively. The pressure amplitudes of 
the LMPE and the HMPE are both 0.1 MPa, and the phase difference 
0 between the pressure oscillations is set at zero. As shown, when 
the frequencies are not of integer multiple relationships, the pres¬ 
sure lift reaches up to 0.18 MPa and is roughly independent of the 
working frequencies. The pressure lifts are nearly zero when the 
frequencies of the LMPE and the HMPE are approximately the same. 
When the frequencies are approximately of other integer multiple 
relationships, the pressure lifts are relatively small. The areas of the 
relatively small pressure lifts around the lines of the integer mul¬ 
tiple relationships decrease with the increase of the multiples, 
which shows that the negative effect of integer multiple relation¬ 
ship of frequency on the pressure lift weakens when the frequency 
difference gets larger. Besides, the negative effects also occur when 
the frequencies are of 5/4,4/3,3/2,5/3, and 5/2 relationships. This is 
because the pressure waveforms of the LMPE and the HMPE have 
less superposition area, and the effective pressure difference and 
opening time of check valve are very limited with these frequency 
relationships when the phase difference is zero. Thus, the working 
frequencies of the engines should avoid the above relationships in 
practice. 
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Fig. 2. Effect of frequency on pressure lift of a CTCS. 
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2.1.2. Effect of pressure amplitude 

Fig. 3 shows the effect of the pressure amplitudes on the pres¬ 
sure lift Ap m for a CTCS. The phase difference 0 is set as zero, and the 
frequencies of the LMPE and the HMPE are 60 Hz and 40 Hz, 
respectively. As shown, the pressure lift is proportional to the 
pressure amplitudes. When the pressure amplitudes are both 
increased from 0.10 MPa to 0.20 MPa, the pressure lift rises from 
0.18 MPa to 0.37 MPa. Moreover, the surface map of the pressure lift 
is distributed symmetrically about the plane of p L ,i = Ph+ which 
indicates that the pressure amplitudes of the LMPE and the HMPE 
have the same effect on the pressure lift. Obviously, more gas is 
pumped from the LMPE to the HMPE when pressure amplitude is 
higher with the present calculation assumptions, which is not 
directly shown in Fig. 3. Thus, the pressure amplitudes should be as 
large as possible in practical applications and methods for 
increasing pressure amplitudes such as adopting pressure amplifier 
[12,28,29] could be adopted. 

Fig. 4 shows the pumping time needed to get the gas pumped 
from the LMPE to the HMPE at different pressure amplitudes. It is 
shown that the pumping time needed is all within several seconds, 
and the maximum average pumping rate of the thermoacoustic 
compression system is up to 5.4 Nm 3 /h. 

2.1.3. Effect of phase difference 

Fig. 5 shows the effect of the phase difference 0 between the 
pressure waveforms on the pressure lift for the CTCS. The pressure 
amplitudes p\ of the LMPE and the HMPE are assumed to be the 
same, and the frequencies are set at 60 Hz and 40 Hz, respectively. 
As shown, the pressure lift reaches the maximum when the phase 
difference 6 is 7t/2, 7tt/ 6 or Htt/6, while the minimum when 0 is 7t/6, 
5tu/ 6 or 3 tu/ 2. For example, when the pressure amplitudes are both 
0.156 MPa, the maximal pressure lift is 1.23 times the minimal one. 
A correlation for the pressure lift could be deduced from Fig. 5 as 

A p m = (|0.3560- sin(1.57r-6»-0.25^)1 + 1.6455)^! -0.01. 

( 8 ) 

The above correlation indicates that the relationship between 
the pressure lift and the phase difference is a full-wave rectification 
waveform with a period of 2 tt/ 3. This is because the transient 
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Fig. 4. Effect of pressure amplitudes on pumping time for a CTCS. 

pressure difference Ap is a function of phase difference with a 
period of 2tt/3. The pressure lift reaches the maximum when the 
peak value of the transient pressure difference Ap is the largest 
with 6 = tu/2 + 2/<tu/3, while the minimum when that of A p is the 
smallest with 0 = tt/6 + 2/cn:/3, where k is an arbitrary integer. The 
phase relationship between the pressure waveforms of the LMPE 
and the HMPE are usually random, so methods to adjust the phase 
relationship between separately running thermoacoustic engines 
should be studied further to realize an efficient phase coupling [30]. 
Besides, the pressure lift is linearly proportional to the pressure 
amplitude for a given phase difference according to the correlation, 
which agrees with the analysis of Fig. 3. Furthermore, the pressure 
lift is found to have a pressure loss of 0.01 MPa in the correlation, 
which is exactly the opening pressure difference Ap c of check valve. 



Fig. 3. Effect of pressure amplitudes on pressure lift for a CTCS. 


Fig. 5. Effect of phase difference on pressure lift for a CTCS. 
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Apparently, the opening pressure difference of check valve should 
be decreased to increase the pressure lift. 

2.14. Transient process simulation of CTCS 

Based on the computation model given before, the dynamic 
thermoacoustic compression process of the CTCS is simulated. For 
the convenience of experimental validation, the geometric and 
operating parameters are set the same as those in the experiments, 
which will be presented systematically later. The experimental 
CTCS is composed of two standing-wave thermoacoustic engines 
and two reservoirs, as shown in Fig. 6. The dimensions of the en¬ 
gines are listed in Table 1. The volumes of LMPR and HMPR are both 
10 L. The calculated working frequencies of the two engines are 
around 65 Hz and 40 Hz when nitrogen is used as the working gas. 

Fig. 7 shows the calculated variation trends of mean pressures 
and instantaneous pumping rates during the compression process 
at different initial mean pressures. As shown, the pressures in 
different parts of the system depart from each other rapidly. It takes 
about 30 s to complete the compression process in the calculation, 
and less than 8 s to complete 90% of the final pressure lift. The 
highest instantaneous pumping rate from the LMPR to the LMPE 
can reach up to 115 Nm 3 /h at 1.80 MPa and 150 Nm 3 /h at 2.40 MPa. 
The calculated mean pressure lift Ap m increases with the initial 
mean pressure, and achieves 0.66 MPa when the initial mean 
pressure is 2.40 MPa. 

2.2. OTCS (open thermo acoustic compression system) 

For a CTCS, gas is always restricted inside it, and gas compres¬ 
sion stops once the pressure difference becomes not high enough to 
open the check valve. Different from the operation of CTCS, OTCS 
continuously pumps gas from a low pressure gas source to a high 
pressure one. The mean pressure in every stage of an OTCS stays 
constant when it is operated steadily, which maintains an effective 
pressure difference to make the compression system work 
continuously. Continuous gas pumping and even novel energy 
conversion systems may be developed based on the principle of 
OTCS. For a steadily operated OTCS, pumping rate is of great 
importance to evaluate the performance, which is quite different 
from the CTCS. Effects of frequency, pressure amplitude, phase 
difference, and mean pressure difference on the average pumping 
rate are calculated and analyzed. 

2.2.1. Effect of frequency 

Fig. 8 shows the effect of the operating frequency on the average 
pumping rate U for an OTCS. The mean pressures pi iYn and pn.m of 


Table 1 

Geometric parameters of the LMPE and HMPE. (D x L: diameter x Length; cp: 
porosity; rh: hydraulic radius). The units for length, diameter and radius are mm. 


Components 

LMPE 



HMPE 



D x L 

(P 

rh 

D x L 

(P 

rh 

Hot duct 

50 x 70 

1 

D/4 

59.5 x 60 

1 

D/4 

Heater 

50 x 80 

0.3 

2 x 10” 3 

59.5 x 50 

0.35 

4 x 10~ 3 

Stack 

50 x 60 

0.84 

3.6 x 10~ 4 

59.5 x 110 

0.4 

3 x 10~ 4 

CHX 

50 x 60 

0.18 

9.5 x 10~ 4 

59.5 x 40 

0.25 

1.5 x 10' 3 

Resonator 

32 x 1050 

1 

D/4 

39 x 2600 

1 

D/4 

Compliance 

81 x 850 

1 

D/4 

85 x 500 

1 

D/4 

cavity 








the LMPE and the HMPE are 2.3 MPa and 2.4 MPa, respectively. The 
pressure amplitudes are both 0.2 MPa and the phase difference 0 
between the pressure waveforms is zero. As shown, the pumping 
rate U is around 50.0 Nm 3 /h at most frequencies, which indicates a 
weak effect of the frequencies on the pumping rate. However, when 
the frequencies are the same, the pumping rate drops down to zero, 
which indicates that no effective superposition area is generated by 
the pressure waveforms at this operating point. The pumping rate is 
low when the frequencies are of 2 and 3 times relationships. For 
example, the pumping rate is about 47.9 Nm 3 /h when the rela¬ 
tionship is 2 times, and it is only 40 Nm 3 /h when the relationship is 
3 times. This is due to the reduced superposition areas of the 
pressure waveforms of the LMPE and HMPE. Different from CTCS, 
when the operating frequency is of larger integer multiple re¬ 
lationships than 3, this influence becomes very weak and even 
vanishes because the superposition area resumes. So, the selection 
of operating frequencies of the thermoacoustic engines should 
avoid the above three integer multiple relationships when 
designing an OTCS. 

2.2.2. Effect of pressure amplitude 

Fig. 9 shows the effect of pressure amplitude on the average 
pumping rate U. The pressure amplitudes of the LMPE and the 
HMPE are assumed to be the same for simplicity. It is clearly shown 
that the pumping rate U rises nearly linearly with the pressure 
amplitude. When the pressure amplitude is increased from 0.1 MPa 
to 0.2 MPa, the pumping rate rises from 16.15 Nm 3 /h to 49.74 Nm 3 / 
h, corresponding to an increase of 208.0%. Therefore, the output 
capability of thermoacoustic engines is critical to the performance 
of the OTCS. Generally speaking, the pressure amplitude of an 
efficient traveling wave thermoacoustic engine is frequently above 
0.3 MPa [15,31,32], which is big enough to generate a pumping rate 
of 72.5 Nm 3 /h. 



Fig. 6. Schematic of experimental setup. P1-P8 denote pressure sensor locations. T1 and T2 denote temperature sensor locations. Gas flows in the direction as shown by arrows. 
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Fig. 7. Variations of mean pressures p m and instantaneous pumping rates U during 
compression at initial mean pressures of (a) 1.80 MPa, (b) 2.13 MPa, and (c) 2.40 MPa in 
calculations. The pressure amplitudes used in calculations are set the same as in the 
experiments before compression process starts. 


2.2.3. Effect of phase difference 

The effect of the phase difference between pressure waveforms 
on the average pumping rate is shown in Fig. 10. The frequencies of 
the LMPE and the HMPE are 60 Hz and 40 Hz, respectively. As 
shown, the variation trend of the pumping rate repeats with a 
period of 2 tt/ 3. When the pressure amplitude is smaller than about 
0.09 MPa, the pumping rate reaches the maximum at a phase dif¬ 
ference of 7 t/2, 7tt/ 6, and Htt/ 6, while the minimum when it is 7 t/6, 
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Fig. 9. Effect of pressure amplitude on average pumping rate for an OTCS. 


5tt/ 6, and 3 tt/ 2. In contrast, the variation trend changes completely 
when the pressure amplitude exceeds 0.09 MPa. Calculation of the 
superposition area at different pressure amplitudes, as illustrated 
by the shadows in Fig. 1, shows that it has the above mentioned 
variation tendency, and further determines the distribution of 
pumping rate in Fig. 10. 

2.2.4. Effect of mean pressure difference 

Fig. 11 shows the effect of the mean pressure difference 
Ap m = p H , m - PL,m on the pumping rate of the OTCS, which can be 
considered as the performance curve like that of traditional pumps 
or fans. The mean pressure of the high-pressure side pn.m is 
assumed to be fixed at 2.4 MPa. As shown, for given pressure am¬ 
plitudes, the pumping rates increases as the mean pressure differ¬ 
ence Ap m decreases. Gas pumping can be easily realized at very low 
pressure amplitude when the mean pressure difference Ap m is 
relatively small. For example, when the mean pressure difference 
Ap m is 0.05 MPa and the pressure amplitude p\ is 0.2 MPa, the 
pumping rate reaches up to 65.37 Nm 3 /h. As the mean pressure 
difference Ap m increases, the pumping rate decreases, and larger 
pressure amplitudes are needed to get the gas pumped. Therefore, 
in order to achieve a certain pressure lift with a relatively high 
pumping rate, a trade-off should be made between them. 



0.050 0 


Fig. 8. Effect of operating frequency on pumping rate for an OTCS. 


Fig. 10. Effect of phase difference on pumping rate for an OTCS. 
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Fig. 11. Effect of mean pressure difference on pumping rate for an OTCS. 

2.2.5. Efficiency analysis of OTCS 

For efficiency analysis, take a specific OTCS as an example, which 
is composed of the two standing-wave thermoacoustic engines 
used in the following experimental study. Fig. 12 shows the effect of 
pressure amplitude on compression efficiency r). The mean pres¬ 
sures of the LMPE and F1MPE are set at 2.3 MPa and 2.4 MPa, 
respectively. The pressure amplitudes of the LMPE and HMPE are 
the same. The phase difference between the engines are set as zero. 
As shown, the compression efficiency reaches the maximum value 
of about 0.57 when the pressure amplitude is 0.092 MPa. When the 
pressure amplitude increases, the difference between the acoustic 
power consumption of the OTCS and the ideal compression work 
becomes larger, and the compression efficiency thus decreases. For 
example, when pressure amplitude is 0.20 MPa, the ideal 
compression work is 46.3 W while that actually required by the 
OTCS is 149.3 W, and the compression efficiency is 0.31. It should be 
noted that the work required by the OTCS is strongly related to the 
thermoacoustic engines selection. If more compact thermoacoustic 
engines utilizing mechanical resonators are used instead of the 
present engines, the actually required acoustic work would be 
remarkably reduced as less acoustic power dissipations occur in¬ 
side the engines. 

Fig. 13 shows the effect of the mean pressure difference between 
the LMPE and the HMPE on the compression efficiency. When the 



P x (MPa) 


Fig. 12. Effect of pressure amplitude on compression efficiency rj. 
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Fig. 13. Effect of mean pressure difference on compression efficiency 77 . 

pressure amplitude is fixed, the compression efficiency first in¬ 
creases and then decreases with the mean pressure difference Ap m . 
When the mean pressure difference is zero, which means the gas is 
not “compressed” from a low mean pressure to a high one, the 
compression efficiency is zero according to its definition. When the 
mean pressure difference is so large that it exceeds the sum of the 
pressure amplitudes of the LMPE and HMPE, for instance, 0.25 MPa 
when pressure amplitude is 0.1 MPa, the efficiency is also zero 
because the pressure of the LMPE can never get higher than that of 
the HMPE and the gas cannot be pumped from the LMPE to the 
HMPE. When the mean pressure difference is fixed, there also exists 
optimal pressure amplitude for the compression efficiency, which 
coincides with the conclusions of Fig. 12. As shown by the dark red 
(in the web version) region in Fig. 13, it is beneficial for the 
compression efficiency when the mean pressure difference and the 
pressure amplitude are both small. The highest compression effi¬ 
ciency reaches about 0.63. 

3. Experimental results 

As shown in Fig. 6, a demonstrative CTCS composed of two 
standing-wave thermoacoustic engines and two reservoirs was 
constructed and tested. The two thermoacoustic engines whose 
working frequencies are about 64 Hz and 40 Hz with nitrogen as 
the working gas are chosen as the LMPE and the HMPE, respec¬ 
tively. The distribution of pressure lift in Fig. 2 shows that the 
pressure lift at these two frequencies is not negatively affected. 
Check valves VI and V3 with an opening pressure of 0.01 MPa were 
supplied by HEROSE (type: 05011.X.0300), and check valve V2 with 
an opening pressure of 0.003-0.021 MPa was supplied by Swagelok 
(type: SS-6C-MM-1/3). Eight calibrated piezoresistive pressure 
sensors with an accuracy of ±0.3% and two calibrated NiCr-NiSi 
thermocouples were installed to measure pressures and heating 
temperatures, respectively. The ball valves were opened one by one 
when the thermoacoustic engines worked steadily. Because of the 
time required to open the ball valves, time lags between them 
occurred during the experiments. 

According to the theoretical analysis, larger pressure amplitudes 
at the two ends of check valves are beneficial to thermoacoustic 
compression and pumping effect. In the experiments, three mea¬ 
sures are taken to obtain as high pressure amplitude as possible. 
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Firstly, the LMPE and the HMPE both work at the maximum 
allowable heating temperature of approximately 650 °C to achieve 
a large temperature ratio between the stack ends. Secondly, con¬ 
necting points for the connecting tubes are set at positions where 
the pressure amplitudes are relatively large. As a result, the con¬ 
necting points for the transfer tube between the engines are set 
near the CHXs, while the ones for the tubes connecting the reser¬ 
voirs are set between the midpoints of the resonators and the CHXs. 
Thirdly, pressure amplifier is adopted to amplify the pressure 
amplitude. The length of the tube connecting the HMPE and the 
check valve V2 is set as 3 m by experimental optimization, which is 
an appropriate length to amplify the pressure amplitude in the 
HMPE without degrading the performance of the thermoacoustic 
engine much. 

Fig. 14(a)-(c) show the variations of mean pressures p m and 
pumping rates U from the LMPR to the LMPE during the ther¬ 
moacoustic compression at different initial mean pressures. The 
pumping rate U is calculated based on the mean pressure variation 
in the LMPR according to the ideal gas equation of state. As shown, 
the mean pressures in the four components depart from each other 
immediately after the compressions start, whose variation trends 




Time (s) 



Time (s) 


Fig. 14. Variations of mean pressures p m and pumping rates U during compression at 
initial mean pressures of (a) 1.80 MPa, (b) 2.13 MPa, and (c) 2.40 MPa in experiments. 
When the initial mean pressure is 2.13 MPa, the opening of check valve V2 is 7 s earlier 
than that of check valve VI, which results in the decrease of the mean pressure of 
LMPE. Thus, larger mean pressure difference between LMPR and LMPE occurs before 
the gas pumping out of LMPR starts. When the initial mean pressure is 2.40 MPa, the 
time gap between the opening of the two valves is only 6 s. Gas pumped out of LMPE is 
not as much as that in 2.13 MPa. As a result, the mean pressure difference between 
LMPR and LMPE is smaller, and the maximum pumping rate in (c) is a little lower than 
that in (b). 


agree well with the transient process simulation in Fig. 7. As stated 
above, the check valves don't start to work simultaneously in 
experimental operations, but one by one from V2, V3 to VI within 
several seconds. As shown in Fig. 14(a), when the check valve V2 is 
first opened, the mean pressure in the LMPE goes down sharply 
while that in the HMPE rises subsequently. Four seconds later when 
the check valve V3 is opened, the gas in the HMPE is pumped into 
the HMPR, resulting in a little decrease of the mean pressure in the 
HMPE and continuous increase of the mean pressure in the HMPR. 
Check valve VI is opened another 4 s later, and the gas in the LMPR 
starts to be pumped into the LMPE. Thus, the mean pressure in the 
LMPE increases a little after that. The change rates of the mean 
pressures and the pumping rate U slow down as the pressure dif¬ 
ferences increase since the superposition areas of the pressure 
waveforms decrease. The variation trends of the mean pressures 
shown in the other two figures are very similar. The average 
pumping rate U reaches up to 4.55 Nm 3 /h during the first 2 s when 
the initial mean pressure is 2.13 MPa, as shown in Fig. 14(b). As the 
time gap between the opening of valve V2 and VI for 2.13 MPa is 1 s 
longer than that for 2.4 MPa, more gas is pumped out of LMPE for 
2.13 MPa, which results in a larger mean pressure difference be¬ 
tween LMPR and LMPE before the gas pumping out of LMPR starts. 
Therefore, the maximum pumping rate for 2.4 MPa is a little lower 
than that for 2.13 MPa. The pressure lift reaches 0.4 MPa when the 
initial mean pressures are 2.40 MPa, as shown in Fig. 14(c). For the 
three cases, the time needed to complete 90% of the final pressure 
lift between the LMPR and the HMPR is less than 125 s. 

The measured and calculated final mean pressure difference 
generated between the LMPR and the HMPR, and the average 
pumping rates during the 90% of the compression processes of the 
above three cases are given in Fig. 15 with respect to the initial 
mean pressure. As shown, the measured pressure amplitudes pn,i 
and pl,i of the HMPE and the LMPE increases with the mean 
pressure. As a result, the final mean pressure lift and average 
pumping rate in the calculation and the experiments both increase 
with the mean pressure. The variation trends of the calculated 
pressure lift and average pumping rate agree well with the exper¬ 
imental ones. When the mean pressure is 2.40 MPa, the pressure 
lifts in the calculation and the experiment are 0.66 MPa and 
0.40 MPa, respectively. The deviations of the average pumping rate 
are much more noticeable. The calculated average pumping rate is 
14.12 Nm 3 /h while the measured one is only 0.51 Nm 3 /h. The de¬ 
viations are mainly due to the pressure waveform deformations 
occurred in experiments, which will be discussed later. As the 
pressure waveforms are so dramatically deformed from sinusoidal 
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Fig. 15. Final mean pressure difference between LMPR and HMPR, and average 
pumping rate during 90% of compression process vs. initial mean pressures. 
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Table 2 

Heating powers of LMPE and HMPE before and after thermoacoustic compression 
starts. 


Initial mean 
pressure 

Heating power before 
compression starts, W 

Heating power after 
compression starts, W 

LMPE 

HMPE 

LMPE 

HMPE 

1.80 MPa 

1016 

845 

909.7 

600.2 

2.13 MPa 

1024 

901.7 

1059 

660 

2.40 MPa 

1092 

912 

1064 

632 


shapes that the effective pressure differences needed to pump gas 
are severely affected. Secondly, the pressure amplitudes in the 
engines are also decreased because of the compression process. 
Finally, the check valves in the real system are different from the 
idealized ones used in modeling. 

The input heating powers of the LMPE and HMPE before and 
after the thermoacoustic compression starts in the experiments are 
given in Table 2. The heating powers are adjusted to keep the 
heating temperatures stable at around 650 °C. As shown, the total 
heating powers are roughly 1.8-2.0 kW before the compression 
starts. After the compression process starts, the heating power 
decreases to about 1.5—1.7 kW, due to the decreases of the pressure 
amplitudes inside the engines. It should be noted that the heating 
power input of the demonstrative system is still large at present, 
due to the poor heat transfer performances of the heater, the 
limited energy conversion efficiencies of the thermoacoustic en¬ 
gines. Besides, the volumes of thermoacoustic engines are rela¬ 
tively large, which results in a large amount of acoustic power 
dissipation. The efficiency of the whole system can be increased if 
much more compact thermoacoustic engines which can provide 
the same pressure amplitude are used. 

Fig. 16 shows the pressure waveforms in the HMPE (PI) and the 
LMPE (P2) when the initial mean pressure is at 2.4 MPa, and the 
timings are denoted accordingly in Fig. 14(c). The gas pumping 
starts at 165 s when check valve V2 is first opened. The pressure 
waveforms are both in good sinusoidal shapes around the same 
mean pressure at 164 s. Once the gas pumping starts, the waveform 
in the HMPE is drastically deformed by the large direct gas flow in 
the first several seconds, as shown at 166 s and 169 s. A little 
deformation of the lower part of the waveform still exists 217 s 
later. When the pressure of the HMPE is higher than that of the 


LMPE, the check valve has no chance to be opened. Thus, the 
deformation of the higher part of the waveform in the HMPE is less 
serious compared to that of the lower part. Almost no deformation 
of the pressure waveform in the LMPE occurs during the whole gas 
pumping process. The reason may be that the check valve V2 is 
much closer to the LMPE and the tube connecting to the HMPE is 
much longer than that to the LMPE. As a result, the long tube 
connected to the HMPE acts as a pressure amplifier and has a strong 
influence on the performance of the HMPE. 

Figs. 17 and 18 show the pressure waveforms before and after 
the check valves VI and V2. As shown in Fig. 17, the pressure P6 in 
the LMPR oscillates at the same frequency as the LMPE after the 
check valve VI is opened, and the phase difference between P8 and 
P6 increases with the time. The mean pressure in the LMPR changes 
very little from 337 s to 397 s, and the pressure oscillation in the 
LMPR almost vanishes at 397 s, which means the gas pumping from 
the LMPR to the LMPE nearly stops even though remarkable pres¬ 
sure difference still exists. A maximum instantaneous pressure 
difference of about 0.07 MPa between P8 and P6 at 397 s indicates 
that the actual opening pressure difference of the check valve VI is 
much larger than the nominal one, which affects the pressure lift 
significantly. The pressure amplitude of P8 rarely has any change 
during the process, indicating the steady performance of the LMPE 
during the gas pumping process. 

Fig. 18 shows the pressure waveforms before (P4) and after (P3) 
the check valve V2. The pressure waveform of P3 after the check 
valve is also seriously deformed as that in the HMPE (PI) as shown 
in Fig. 16. Besides, the pressure waveform of P3 is more seriously 
influenced at the lower part when the gas is sucked from the LMPE 
side into the HMPE side. Moreover, the peak-to-peak amplitude of 
P3 is about 0.31 MPa while that of PI is only 0.18 MPa at 217 s, 
which indicates the pressure amplification effect of the pressure 
amplifier. The pressure waveform of P4 is also slightly deformed, 
due to the much higher direct gas flow rate within the thin con¬ 
necting tube than that in the LMPE. In brief, the gas pumping re¬ 
sults in the pressure waveform deformations in the thermoacoustic 
compression system, which is not taken into account in previous 
calculation yet. The deformed pressure waveforms then decrease 
the effective pressure differences across the check valves, and affect 
the gas pumping rate negatively in turn. Thus, these two factors 
lead to the formation of a negative feedback system, and make the 
gas pumping rate fall in a reasonable range. 
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Fig. 16. Pressure waveforms in HMPE and LMPE during thermoacoustic compression process. The gas pumping starts at 165 s when check valve V2 is opened. 
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Fig. 17. Pressure waveforms before and after check valve VI during thermoacoustic compression process. Check valve VI opens at 172 s. 
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Fig. 18. Pressure waveforms before and after check valve V2 during thermoacoustic compression process. Check valve V2 opens at 165 s. 
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Till now, the key factors that affect the performance are 
revealed. It is shown that the performance can be further improved 
in the following aspects. Firstly, more powerful, efficient and 
compact thermoacoustic engines should be designed to enlarge 
pressure amplitudes and decrease acoustic power losses. Pressure 
amplifier can be further optimized to amplify pressure amplitudes 
without affecting the thermoacoustic conversion too much. Sec¬ 
ondly, the flow resistance, the opening pressure difference, and the 
response time of the check valves should be minimized to reduce 
undesirable losses. Thirdly, effective approaches to weaken the 
influence of pressure waveform deformations before and after 
check valves need to be explored. Last but not least, the coupling 
mechanism of thermoacoustic engines and check valves should be 
further investigated. Efficient thermoacoustic conversion exposed 
both to alternating and direct gas flows should be studied. 

4. Conclusions 

We have proposed a thermoacoustic compression effect based 
on the pressure oscillation nature of thermoacoustic engines and 


the gas flow rectification effect of the check valves. On the basis of 
the effect, closed and open thermoacoustic compression systems 
can be realized to get gas compressed and pumped. Theoretical 
analyses and computational results have shown the application 
potential of such systems for its remarkable merits, such as large 
pumping rate, high pressure lift, simple configuration, thermally 
driven nature, and oil-free operation. 

Computation shows that the operating frequencies of the ther¬ 
moacoustic engines should not be of small integer multiple re¬ 
lationships for efficient operation. Large pressure amplitudes are 
beneficial to the performance of both the closed and the open 
systems. It is also shown that pressure amplitude and mean pres¬ 
sure difference between the low and high mean pressure engines 
should be optimized to increase the compression and pumping 
efficiency. Besides, it is beneficial for the compression efficiency if 
more compact thermoacoustic engines with high pressure ampli¬ 
tudes are used in an open system. A demonstrative closed ther¬ 
moacoustic compression system are then experimentally studied. 
The variation trends of the mean pressures during the compression 
process, final pressure lifts, and average pumping rate in the 
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calculations and the experiments agree well. In the experiments, a 
maximum average pumping rate of 4.55 Nm 3 /h and a maximum 
pressure lift of 0.4 MPa were achieved when the initial mean 
pressures were 2.13 MPa and 2.40 MPa, respectively. Furthermore, 
it is found that the direct gas flow can cause serious pressure 
waveform deformation, which then affects the gas pumping 
negatively in turn. 

Acknowledgments 

This research is financially supported by the Major State Basic 
Research Development Program of China (973 Program) under 
contract No. 2010CB227303 and National Natural Science Founda¬ 
tion of China under contract No. 61077035. 

References 

[1 ] McMahon HO, Gifford WE. A new low-temperature gas expansion cycle-part I, 
in advances in cryogenic engineering. New York: Plenum Press; 1960. 
p. 354-66. 

[2] Sun DM, Xu Y, Chen HJ, Shen Q Zhang XJ, Qiu LM. Acoustic characteristics of a 
mean flow acoustic engine capable of wind energy harvesting: effect of 
resonator tube length. Energy 2013;55:361-8. 

[3] Sun DM, Xu Y, Chen HJ, Wu K, Liu KK, Yu Y. A mean flow acoustic engine 
capable of wind energy harvesting. Energy Convers Manag 2012;63(0):101—5. 

[4] Yu Y, Sun DM, Wu K, Xu Y, Chen HJ, Zhang XJ, et al. CFD study on mean flow 
engine for wind power exploitation. Energy Convers Manag 2011;52(6): 
2355-9. 

[5] Slaton WV, Zeegers JCH. An aeroacoustically driven thermoacoustic heat 
pump. J Acoust Soc Am 2005;117(6):3628—35. 

[6] Slaton WV, Zeegers JCH. Acoustic power measurements of a damped aero¬ 
acoustically driven resonator. J Acoust Soc Am 2005;118(1):83—91. 

[7] Ullmann A. The piezoelectric valve-less pump—performance enhancement 
analysis. Sens Actuators A Phys 1998;69(1):97—105. 

[8] Ullmann A, Fono I. The piezoelectric valve-less pump-improved dynamic 
model. J Microelectromech Syst 2002;11(6):655—64. 

[9] Lee DG, Or SW, Carman GP. Design of a piezoelectric-hydraulic pump with 
active valves. J Intell Mater Syst Struct 2004;15(2):107—15. 

[10] Swift GW, Backhaus S. A resonant, self-pumped, circulating thermoacoustic 
heat exchanger. J Acoust Soc Am 2004;116(5):2923—38. 

[11] Gao B, Wu ZH, Luo EC, Dai W. Experimental demonstration of a novel heat 
exchange loop used for oscillating flow systems. 985. Advances in cryogenic 
engineering, vol. 53a and 53b; 2008. p. 1114—21. 

[12] Sun DM, Wang K, Xu Y, Shen Q Zhang XJ, Qiu LM. Thermoacoustic 
compression based on alternating to direct gas flow conversion. J Appl Phys 
2012;111(9). 


[13] Sun DM, Wang K, Qiu LM. Multistage thermoacoustic compression system. 
Chinese Patent, Patent No. 200810163521.3. 

[14] Sun DM, Wang K, Lou P, Zhao YT, Qiu LM, A refrigeration system driven by 
thermoacoustic compressor. Chinese Patent, Application No. 201110261644.2. 

[15] Backhaus S, Swift GW. A thermoacoustic Stirling heat engine. Nature 
1999;399(6734):335—8. 

[16] Bisio G, Rubatto G. Sondhauss and Rijke oscillations — thermodynamic anal¬ 
ysis, possible applications and analogies. Energy 1999;24(2):117—31. 

[17] Chun W, Oh SJ, Lee YJ, Lim SH, Surathu R, Chen K. Acoustic waves generated 
by a TA (ThermoAcoustic) laser pair. Energy 2012;45(1):541—5. 

[18] Zhao D, Ji C, Li S, Li J. Thermodynamic measurement and analysis of dual¬ 
temperature thermoacoustic oscillations for energy harvesting application. 
Energy 2014;65(0):517-26. 

[19] Flitcroft M, Symko OG. Ultrasonic thermoacoustic energy converter. Ultra¬ 
sonics 2013;53(3):672—6. 

[20] Roberts DC, Li H, Steyn JL, Yaglioglu O, Spearing SM, Schmidt MA, et al. 
A piezoelectric microvalve for compact high-frequency, high-differential 
pressure hydraulic micropumping systems. J Microelectromech Syst 
2003;12(1):81—92. 

[21] Park JM, Taylor RP, Evans AT, Brosten TR, Nellis GF, Klein SA, et al. 
A piezoelectric microvalve for cryogenic applications. J Micromech Microeng 
2008;18(1):015023. 

[22] Gradin H, Clausi D, Braun S, Stemme G, Peirs J, van der Wijngaart W, et al. 
A low-power high-flow shape memory alloy wire gas microvalve. 
J Micromech Microeng 2012;22(7):075002. 

[23] Gedeon D. Sage stirling-cycle model class reference guide. 5th ed. Athens, OH: 
Gedeon Associates; 2007. 

[24] Wollan JJ, Swift GW, Backhaus SN, Gardner DL. Development of a thermoa¬ 
coustic natural gas liquefier. In: AIChE Meeting; March 11—14, 2002 [New 
Orleans LA]. 

[25] Swift GW. Thermoacoustics: a unifying perspective for some engines and 
refrigerators. Sewickley, PA: Acoustical Society of America; 2002. 

[26] Qiu LM, Lai BH, Zhao YT, Sun DM, Zhang XJ, Li YF. Study on the onset tem¬ 
perature of a standing-wave thermoacoustic engine based on circuit network 
theory. Sci China Technol Sci 2012;55(10):2864—8. 

[27] Lai BH, Qiu LM, Li YF, Lou P, Sun DM. Simulation of onset process of standing 
wave thermoacoustic engine based on thermoacoustic network theory. 
J Zhejiang Univ Eng Sci 2011;45(6):1130—5 [in Chinese]. 

[28] Sun DM, Qiu LM, Wang B, Xiao Y. Transmission characteristics of acoustic 
amplifier in thermoacoustic engine. Energy Convers Manag 2008;49(5): 
913-8. 

[29] Sun DM, Qiu LM, Wang B, Xiao Y. Novel Helmholtz resonator used to focus 
acoustic energy of thermoacoustic engine. Appl Therm Eng 2009;29(5—6): 
945-9. 

[30] Gillman BJ. A study of coupled thermoacoustic engines. Doctor of Philosophy. 
The University of Utah; 2009. 

[31 ] Sun DM, Qiu LM, Zhang W, Yan WL, Chen GB. Investigation on traveling wave 
thermoacoustic heat engine with high pressure amplitude. Energy Convers 
Manag 2005;46(2):281-91. 

[32] Yu GY, Luo EC, Dai W, Wu ZH. An energy-focused thermoacoustic-Stirling heat 
engine reaching a high pressure ratio above 1.40. Cryogenics 2007;47(2): 
132-4. 


